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OBJECTIVES 

The moving molten contact line plays a key role 
in a number of applications, including certain coating 
processes, spin casting of metals, and various droplet 
based manufacturing techniques. The purpose of our 
program is to develop a basic understanding of how a 
molten material front spreads over a solid substrate 
that is below its melting point, arrests, and freezes. 
We are concerned particularly with the dynamic con- 
tact line that occurs in molten droplet deposition on 
cold substrates. Our hope is that the work will provide 
a scientific knowledge base for certain new applica- 
tions, such as the “printing” of arbitrary three- 
dimensional objects by precise deposition of individual 
molten microdrops (Gao & Sonin, 1994), and at the 
same time contribute generally to the understanding of 
low Bond number, capillarity-driven liquid dynamics 
problems. 

The main obstacle to progress at this time is our 
lack of knowledge of the basic processes in molten 
droplet deposition. We do not yet know the basic laws 
that govern the dynamics of the molten contact line 
and its eventual arrest by freezing, and consequently 
are not in a position to answer the most fundamental 
question in molten droplet deposition: given the melt 
and target materials, how do we choose the droplet 
size and deposition conditions so as to produce a suit- 
able solid shape? The answer to this question depends 
on the motion and arrest of the molten contact line. 


ARREST OF A MOLTEN CONTACT 
LINE MOVING OVER A SUBCOOLED 
SOLID OF ITS OWN KIND: 
EXPERIMENTS AND THEORY 

A molten droplet usually spreads over a cold sub- 
strate much quicker than it solidifies in bulk, that is, 
the drop remains essentially in liquid form until the 
contact line arrests. The point of contact line arrest 
determines the final shape of the droplet. An ordinary 
liquid’s contact line arrests at a certain “equilibrium” 
contact angle. When a molten material is involved, 
however, the spreading is typically arrested by contact 
line freezing. Indeed, a little reflection shows that in 


homologous deposition (melt and target of the same 
material, one above and the other below the freezing 
point), there exists no meaningful final equilibrium 
contact angle. 

Using an experiment where the rate of contact line 
advance was externally controlled, Schiaffino and 
Sonin (1997a) investigated the arrest of the contact 
line of a molten microcrystalline wax spreading over a 
subcooled solid “target” of the same material. They 
found (Fig. 1) that contact line arrest takes place at an 
apparent liquid contact angle which depends primarily 
on the Stefan number S=c(T r T l )/L based on the tem- 
perature difference T f -T t between the fusion point and 
the initial target temperature, c being the specific heat 
and L the latent heat of fusion. Neither the melt’s su- 
perheat nor the prior history of the contact line motion 
appeared to significantly affect the contact angle at the 
point where the contact line comes to a halt. 



Schiaffino and Sonin proposed the following 
physical mechanism for contact line arrest due to 
freezing: arrest occurs when the liquid’s dynamic 
contact angle approaches the angle of attack 0 s of the 
solidification front just behind the contact line (Fig. 2). 
Contact line advance is possible only as long as the 
liquid’s contact angle 0 a exceeds 0 s . 

In a second paper Schiaffino and Sonin (1997b) 
showed that the conventional continuum equations and 
boundary conditions for the slow spreading of a pure 
molten material over a solid of its own kind have no 
meaningful solution for 0 S . The quantity 0 S is deler- 
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mined by the heat flux just behind the contact line, and 
the heat flux in the conventional mathematical model 
is singular at the contact line. However, by comparing 
experimental data with numerical computations, 
Schiaffino and Sonin estimated that for microcrystal- 
line wax, the breakdown of the conventional solidifi- 
cation model occurs within a distance of order 0.1-1 
pm of the contact line. The physical mechanism for 
this breakdown is as yet undetermined, and conse- 
quently no first-principles theory exists at this time for 
e 5 and thus for contact line arTest. However, the func- 
tional relationship 0 s =f(S) can in principle be deter- 
mined experimentally for each material. 



wax contact line reference frame. 


MOLTEN DROPLET DEPOSITION: 
SIMILARITY LAWS AND 
EXPERIMENTS 

A third paper by Schiaffino and Sonin (1997c) 
presents a study of the deposition process of molten 
droplets, from the instant of contact to complete solidi- 
fication. This study lays down a framework for under- 
standing low to moderate Weber number molten 
droplet deposition in terms of similarity laws and ex- 
perimentation. The study is based on experiments 
with three molten materials — molten wax on solid 
wax, water on ice, and mercury on frozen mer- 
cury — which between them span a considerable range 
of the deposition/solidification similarity parameters. 
Based on experiments from the highly viscous limit to 
the inertia dominated limit, correlations are obtained 
for the spreading velocity, spreading time scales, the 
spreading factor (i.e. ratio of deposited drop’s final 
footprint radius and the drop’s initial radius), post- 
spreading liquid oscillation amplitudes and time scales, 
and bulk solidification time scales. At low to moderate 
Weber numbers, spreading is driven principally by 


interfacial forces at the contact line rather than by im- 
pact phenomena. The main dynamic similarity pa- 
rameter is the Ohnesorge number Z=ji/(pGa) , where 
ji is the absolute viscosity, p is the density, a in the 
surface tension, and a is the drop radius. The principal 
similarity parameter for contact line arrest is the Stefan 
number. 

FURTHER EXPERIMENTS 

Experiments have been carried out to (i) extend 
the investigation of contact line freezing to different 
materials and conditions, (ii) investigate the relation- 
ship between the dynamic contact angle and the con- 
tact line speed during spreading, and (iii) better define 
the thermal and solidification properties of the micro- 
crystalline wax we have been using. 

Contact line motion and arrest were studied by 
depositing molten drops with radius of order 1 mm on 
solid targets below the melt’s fusion point, and photo- 
graphing the process using high-speed video micros- 
copy. Complete spreading, arrest and solidification 
events were recorded for molten wax on solid wax, for 
water on ice, and for molten wax deposited on inclined 
wax surfaces in order to investigate gravitational ef- 
fects on molten contact line dynamics and arrest (Fig. 
3). From the high-speed video data we obtain empiri- 
cal relationships between the instantaneous values of 
the capillary number Ca= (lU/o (U is the contact line 
advance speed) and melt contact angle 0 a . This rela- 
tionship is discussed in the next section. 

Post-arrest angles were obtained for molten solder 
on glass, molten solder on solid solder, and for octaco- 
sane on solid octacosanc with two different (measured) 
surface roughnesses. Octacosanc (C 28 H 5g ) is a low 
melting point alkane with a distinct fusion point and 
latent heat release, and has well known thermal prop- 
erties. We also investigated molten octacosane on 
clean glass and on glass coated with two types of self- 
assembled monolayers, one of which chemically re- 
sembles octacosane. These results will be discussed at 
the conference. 

In addition, a study was completed of the thermal 
and solidification properties of microcrystalline wax, 
including its thermal conductivity k and enthalpy h(T) 
and specific heat c(T) at temperatures between 0°C and 
115°C (Torrcsola, 1998). This wax has a distinct fu- 
sion point at 90 C (Schiaffino, 1996), but releases 
latent heat over a broad temperature range (about 
20 C) below that point. 
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THEORY FOR VELOCITY VS. ANGLE 
AT A MOLTEN CONTACT LINE 

A significant body of both experimental and theo- 
retical work has built up over the past two decades on 
the relationship between speed and apparent contact 
angle at an ordinary liquid contact line. See for exam- 
ple the reviews by Dussan V. (1979), dc Gennes 
(1985), Kistler (1993), Blake (1993), and the com- 
ments in Schiaffino & Sonin (1997a). One of the sim- 
plest approximate forms of this relationship is the 
Hoffman-Tanncr-Voinov law for small Ca (Hoffman, 
1975; Tanner, 1978; Voinov, 1976, 1978; see also 
Boendcr, et al , 1991), 

c.= K$-ei) (i) 

where 0 m is the liquid’s contact angle at the molecular 
scale, and is usually identified with the equilibrium 
angle. Hoffman’s experimental value for k is about 
1.3xlO' 3 rad 3 , although it does in fact depend weakly 
on flow system size. 

Other than some results mentioned in Schiaffino 
and Sonin (1997a), no work has previously been done 
on the Ca vs. 0 a relationship for a molten contact line. 
The molten dynamic contact line problem differs sig- 
nificantly from the ordinary one in that a solidification 
front forms underneath the molten material as it moves 
across the subcooled solid, and the material stream- 
lines move into the solid as sketched in Fig. 2. Schiaf- 
fino & Sonin (1997a,b) argued that the solidification 
front at the contact line is wedge-like with a slope 0 s 
relative to the “horizontal”, and deduced 0 s (S) empiri- 
cally. 

We have adapted Voinov’s methodology to the 
molten contact line and obtained solutions for the liq- 
uid flow field near the contact line and for Ca vs. 0 a 
(Duthaler, 1998). Since the actual solidification front 


shape near the contact line is unknown, we approxi- 
mate it as piecewise continuous, with constant slope 0 s 
up to a cut-off distance r=X and zero slope thereafter 
(Fig. 4). The flow in the liquid region is assumed to be 
inertia-free and quasi-steady near the contact line, with 
viscous forces balancing the pressure gradient result- 
ing from interfacial curvature. In general the equa- 
tions require a numerical solution, but in the final 
stages of the arrest process, where (0 a -0 5 )«0 s , they 
yield the following surprisingly simple analytical ex- 
pression 


Ca = -^-^e a - 0 s ) 4 -( e m - 0 s) 4 } (2) 



Figure 4 Piecewise continuous solidification front geometry used in 
modeling the molten dynamic contact line problem 


Here, k plays a role similar to the coefficient in Eq. 
(1), but docs not have exactly the same value. Eq. (1) 
is in essence an analog of the Hoffman-Tanner-Voinov 
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law for molten materials approaching the final stage of 
arrest. Note that near contact line arrest, Ca scales as 
the 4th power of the (small) difference 0 a -0 s . This be- 
havior differs dramatically from that of the Hoffman- 
Tanner-Voinov Law, Eq. (1). 

The theory can be applied, for example, to calcu- 
late the time scale associated with a molten drop 
spreading event. For example, in the limit of small 
Bond and Capillary numbers and small interfacial an- 
gles (0 a ,0 m , 0 S «1), we obtain 

1 = (3) 

( 0 a - 0 s ) 3 _ ( 6 o- 0 s ) 3 Sin 0 s T 


where 0 o is the value of 0 a at t=0, 


x = 


JL 

3kg 



(4) 


Before a useful theory for molten contact line ad- 
vance and freezing can be developed, we need a real- 
istic model for the heat flux at the cont act line. One of 



Figure 5 0 8 vs. Ca relationship for a molten drop of microcrystal- 
line wax spreading over a solid target of the same kind. The numeri- 
cal solution begins when R/a=l .64 and H/a=0.90. 


is a characteristic time, and V denotes the volume of 
the drop. Note that 0 a -0 ~ t' 1 " for this case. Applying a 
similar scaling analysis to Schiaffino’s force feed ex- 
periment (see Schiaffino and Sonin, 1997b) shows that 
very near the arrest of a force fed droplet, 0 a -0 s ~ f ut . 

The magnitude and trends predicted by the new 
spreading model are in basic agreement with our ex- 
perimental results (Fig. 5). The plot shows empirical 
data for a drop deposition characterized by S=0.15, 
along with the predictions of various spreading mod- 
els. The solid line represents the numerical solution of 
the new hydrodynamic equations derived for melts, 
based on applying Voinov’s methodology to the sim- 
plified, wedge-like solidification interface of Fig. 4. 
We have set the wedge size A. equal to the cut-off dis- 
tance computed by Schiaffino and Sonin (1997b), 
k=0.1 pm. We also take 0=20°, and assume that 0 m =0 s 
and that the bulk of the drop remains hemispherical 
during the final, slow stage of spreading. The ap- 
proximate solution given by Eq. (2) with 0 =0^=20* is 
shown using long dashes. Since no a priori model for 
0 s exists, we choose the value 0 s =2O* to ensure good 
agreement with the experimental final arrest angle, 
28.5*. For reference, Fig. 5 also shows Eq. (1) with 0 m 
set equal to the final arrest angle. At Ca<10 2 this 
equation overestimates Ca for given 0 a by nearly an 
order of magnitude. 

HIGH-RESOLUTION MICROSENSORS 
FOR THE THERMAL PROCESSES 
NEAR THE MOVING CONTACT LINE 


the key objectives of our program is the development 
of micron-scale sensors for measuring the transient 
temperature at a point on the substrate surface as a 
molten contact line moves over it. The hope is that 
from this temperature history, we can back out infor- 
mation on the heat flux distribution near the contact 
line and obtain important clues about the nature of the 
thermal process in that region. 

The sensors arc of the thermistor type, microfabri- 
cated with silicon-based technology on either pure 
silicon or amorphous silicon dioxide wafers, 100 mm 
in diameter and 0.5 mm thick. Each wafer has room 
for 60 chips. Each chip has 32 sensors on its surface, 8 
arrayed in a cross-pattern and spaced 400 pm apart, 8 
arrayed in line with 200 pm spacing, and 16 arrayed in 
line with 50 pm spacing (Figs. 6 and 7). The sensors 
in the two coarsely spaced arrays are 2.5 pm square, 
and those in the fine array arc 1.5 pm square. The 
time response is better than 10 ps. The sensing ele- 
ments are made of boron-doped silicon and are electri- 
cally connected to the outside circuitry by aluminum 
circuit traces and gold wires. Temperature is sensed 
via Wheatstone bridges which detect the sensors’ re- 
sistance variations due to temperature. Each individ- 
ual sensor must be separately calibrated. Data is ac- 
quired with a PC 16-channel interface board, at a 
maximum rate of 10 6 samples/sec. 

At the time of writing, the chips have been fabri- 
cated on wafers (we have several hundred chips of 
each of the two substrate materials), and a dozen chips 
have been cut and packaged. The data acquisition 
system has been constructed, and sensor calibration 
and testing is in progress. Some results on thermal 


145 



transients during contact line passage will be available 
for the August conference. While we expect that the 
data will provide information on the near-contaci-linc 
heat transfer process, we also foresee possible prob- 
lems. First, small as they are, the sensors are actually 
slightly larger than the cutoff lengths computed by 
Schiaffino and Sonin (1997b), and the spatial resolu- 
tion may be insufficient for resolving the near-contact- 
line region where the conventional model departs from 
reality. Second, due to the microfabrication method 
the conductors protrude about 0.5 pm above the sub- 
strate surface, which may be large enough to affect 
contact line motion. Finally, a sensor’s temperature 
history will depend not only on the heat flux distribu- 
tion into it from the fusion front, but also on the ther- 
mal properties of the substrate below it and the solidi- 
fied melt between it and the fusion front (Fig. 2). The 
heat flux distribution in the contact line region must 
therefore be unfolded from computations of the overall 
system’s transient response. Whether these potential 
problems arc significant will become clear as the work 
proceeds. 
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Figure 6 Silicon thermal sensor chip showing the sensor arrays 



Figure 7 Close-up of the coarse cross-pattern sensor array 
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